Competition between isoscalar and isovector pairing correlations in N = Z nuclei 
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We study the isoscalar (T = 0) and isovector (T = 1) pairing correlations in N — Z nuclei. They 
are estimated from the double difference of binding energies for odd-odd N — Z nuclei and the 
odd-even mass difference for the neighboring odd-mass nuclei, respectively. The empirical and BCS 
calculations based on a T = and T = 1 pairing model reproduce well the almost degeneracy of 
the lowest T = and T — 1 states over a wide range of even-even and odd-odd N = Z nuclei. It is 
shown that this degeneracy is attributed to competition between the isoscalar and isovector pairing 
correlations in N — Z nuclei. The calculations give an interesting prediction that the odd-odd 
N = Z nucleus 82 Nb has possibly the ground state with T = 0. 
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There is a current topic with increasing interests in 
studying isovector (T = 1) and isoscalar (T = 0) proton- 
neutron (pn) pairing correlations in N — Z nuclei 1]. 
At present, it is not clear whether pn pairing correla- 
tions are strong enough to form a static condensate. It 
is well known that an experimental signature of likc- 
nucleon proton-proton (pp) and neutron-neutron (mi) 
J = pairing correlations in nuclei with neutron excess 
is the odd-even mass difference, which is extra binding 
energy of even-even nuclei relative to that of odd-mass 
nuclei. However, the odd-even mass differences for even- 
even N = Z nuclei are larger than those of the neigh- 
boring even-even N = Z + 2 nuclei, and it reflects the 
gain in pairing due to stronger pn correlations • It has 
recently been shown 0, Q that the three-point odd-even 
mass difference for an odd-mass nucleus with neutron 
excess is an excellent measure of pp and nn pairing cor- 
relations in neighboring even-even nucleus, although it is 
still controversial [jj. This conclusion suggests that the 
pp and nn pairing correlations in N = Z even-even nuclei 
also can be estimated from the odd-even mass difference 
of neighboring odd-mass nuclei with N = Z + 1 . On the 
other hand, the pn pairing can be estimated from the 
double difference of binding energies Q. When we as- 
sume isospin symmetry in N ss Z nuclei, the T = 1 pn 
pairing and like-nucleon [pp and nn) pairing are classified 
in the same T = 1 pairing correlations, and the former 
correlation energy should be the same as the latter one. 

Odd-odd N = Z nuclei are an ideal experimental lab- 
oratory for the study of pn pairing correlations. It is well 
known that the lowest T = and T = 1 states in odd-odd 
N = Z nuclei are almost degenerate and exhibit the in- 
version of the sign of the energy difference Et=i — -Er=o> 
while all even-even N = Z nuclei have the T = ground 
states and the T = 1 excited states with l arge excita- 
tion energies. Several authors 0,0 HUGH E3 already 
pointed out that this degeneracy in odd-odd N = Z nu- 
clei reflects the delicate balance between the symmetry 



energy and the pairing correlations. The T = and 
T = 1 ground-state binding energies of N — Z nuclei 
were calculated by using an algebraic model based on 
IBM-4 [13 In this paper, we study the T — and 
T = 1 pairing correlations from a phenomenological point 
of view, and analyze them in the BCS calculations within 
a schematic model that includes T = 1 and T = pairing 
interactions. 

We begin with the estimation of T = 1 pairing corre- 
lations in N = Z nuclei. A typical indicator for T = 1 
pairing correlations is the following three-point odd-even 
mass difference: 

(— 1 \ N 

A^(Z.N) = ^—^ [B (Z,N + 1) 

- 2B(Z,N) +B(Z,N- 1)], (1) 

where B(Z, N) is the negative binding energy of a sys- 
tem. Since B(Z, N ± 1) w B(Z, N) + A ± A based 
on standard BCS theory with pairing gap A leads to 
An (Z, N) w A, the indicator A^ is often interpreted 
as a measure of the empirical pairing gap. However, 
it is well known that values of [Z = even, N) are 
large for even-A^ and small for odd-TV. It was dis- 
cussed 3] that A„ (Z = even, N = odd) is an excellent 
measure of T = 1 pairing correlations, and the diffcr- 

(3) 

ences of A„ at adjacent even- and odd-N nuclei reflect 
the mean- field contributions. From a view point of the 
semi-empirical mass formula, the above indicator is well 
known to be affected by the symmetry energy term in 
the liquid-drop model. In the macroscopic-microscopic 
shell model, however, the curvature contribution cancels 
out the symmtery energy contribution as pointed out by 
Satulta et al.^. What does the magnitude of the pair- 
ing gap in the N — Z nuclei mean? We suggest that 
A'n'(Z,Z+l) of odd -mass nucleus should be regarded 
as pure pairing gap in N = Z adjacent even-even and 
odd-odd nuclei. For the N = Z nuclei, the four and 
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FIG. 1: The experimental odd-even mass differences 
An (Z = even, Z + 1) (solid diamonds) in odd-mass nuclei 
with N = Z + 1 , and the pairing gaps (open circles) obtained 
by the BCS calculations. The solid curve is 5.18A _1/3 and 
the dashed curve denotes I2A -1 / 2 . 



five point indicators cannot be adopted because they in- 
clude large contributions from mean filed and pn cor- 
relations [3, @ • Figure 1 shows experimental values of 



A^ in odd-mass nuclei, where we plot AX'(Z,Z + i) for 
16 < A < 60. When there is no data of A { n\z, Z + 1) 
for 60 < A < 110, we adopt for nearest nuclei with 
N = Z + 1. The expected quenching of neutron pairing 
at magic (or semi-magic) particle number N or Z=14, 
28, 40, and 50 is clearly seen in the figure. The standard 
curve 12A" 1 / 2 is also shown as a guide eye in Fig. 1. 
We can see that the average pairing gap is smaller than 
the values of the curve 12A~ 1//2 . The global trend can 
be fitted by the curve 5.18A -1 / 3 MeV, as discussed in 
recent analyses where T = 1 pairing gap At=i 

obtained from some binding energy difference is fitted by 
the mass-dependence A -1 / 3 different from the standard 
one YlA~ x l 2 . The difference between the two curves is 
quite large for light nuclei, while it is small for heavy 
nuclei. The average gap was recently analyzed 0] by 
A = a + pA^ 1 / 3 which has theoretical foundation. This 
analysis also supports the weaker mass- dependence. We 
now consider the following pairing Hamiltonian to de- 
scribe the T = 1 pairing correlations: 



H = Hn 



(2) 



where e a is the single-particle energy and P K is the J = 
pair operator with isospin T = 1, T z = k. Implying 
isospin invariance to the above Hamiltonian, the pairing 
part Hp includes the isovector pn interactions. The stan- 



dard BCS calculations with the pairing Hamiltonian (j2J 
were performed in sd and fpg shells. We adopted single- 
particle energies from a spherical Woods-Saxon poten- 
tial in the BCS calculations. The pairing force strength 
G = 24.5 /A was chosen so as to fit the experimental odd- 
even mass difference AiT (Z — even, Z + 1) in odd-mass 
nuclei. The BCS results for A > 40 almost agree with 
the experimental odd-even mass differences, and more- 
over reproduce the shell effects. The BCS calculations 
reproduce well the behavior of the observed odd-even 
mass difference over a wide range of N — Z nuclei. Thus 
the T = 1 pairing correlations can be estimated from 
the odd-even mass difference A„ (Z = even, Z + 1) in 
odd-mass nuclei. 

To describe the pn pairing correlations in odd-odd N — 
Z nuclei, let us estimate the following double difference 
of binding energies 0, El : 

Aj„(Z,iV) = ±[B(Z,N) T -B(Z,N-1) 

- B{Z -1,N)+B{Z -1,N -I)], (3) 

where B(Z,N) T is the binding energy of lowest state 
with isospin T in odd-odd N = Z nuclei. Figure 2 shows 
the double difference of binding energies calculated from 
the experimental binding energies. The odd-even mass 
differences for odd-mass nuclei are also displayed. Then 
we can see that the (Z = even, Z+l) agrees with the 
A?'= 1 (Z + 1,Z+ 1). This means that T = 1 pn pairing 
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FIG. 2: The pn pairing gaps estimated from the double dif- 
ferences of experimental binding energies. The solid circles 
denote the T — pn pairing gap, and the solid triangles the 
the T — 1 pn pairing gap. The odd-even mass differences 
in odd-mass nuclei with N = Z + 1 are shown by the open 
squares. The dashed curve is the half of the T = pairing 
force strength k°. 
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for odd-odd N = Z nuclei have the same correlation 
energy as the like-nucleon nn pairing, A„ = Aj^ 1 , when 
assuming isospin symmetry. Thus, the indicator Aj" 1 
gives the T = 1 pn pairing gap in N = Z nuclei. The 
Apn° can be regarded as the T — pn pairing gap as 
well. Figure 2 with these estimations indicates that the 
T = pn correlations are superior to the T = 1 pn 
correlations in the ground states of sd shell nuclei, and 
the inversion occurs in the pf shell nuclei. The T = 
pn pairing gap Aj^ cannot be explained by the T = 1 
pairing Hamiltonian 

In a previous paper 2], it has been shown that the 
T = matrix elements of the monopole field V^(a,b) 
are significantly larger than the T = 1 ones, and are very 
important in determining the double differences of bind- 
ing energies, where a, b are the single particle orbitals. 
We can see that the matrix elements are quite large for 
isoscalar components but small for isovector components. 
In the USD interaction, the monopole matrix elements 
with T = have values around -3 MeV and are strongly 
attractive. If we assume that the T = monopole matrix 
elements are equal and independent of angular momen- 
tum J and the single particle orbitals, V£~° is reduced 
to the J-independent isoscalar p-n pairing interaction. 
Neglecting T = 1 monopole components, let us add the 
J-independent T = pn pairing interaction 0, Il7| to the 
pairing Hamiltonian 



TT=0 



H + H P -k°J2J2 At JMm(ab)A 



JM.00 



(ab)(A) 



a>b J,M 



where A JM 00 (ab) is the pair operator with spin J and 
isospin T = 0. The T = 1 pairing interaction does not 
contribute to the double difference of binding energies 
Aj=°, and Aj=° « fc°/2. Then, the T — pairing force 
strength fc° = 244.5(1 - 1.67 A" 1 / 3 )/ A is chosen so as 
to fit the T = pn pairing gap as seen in Fig. 2. The 
isovector monopole components in USD are small, except 
for U T j =1 (si/2, S1/2). The deviations from the curve k°/2 
for 30 P and 34 C1 in Fig. 2 would be attributed to the 
large value of isovector component V^ =1 (s 1 / 2 , s 1 / 2 ). We 
recently introduced |l7| monopole corrections to improve 
the energy levels of 48 Ca, etc. In this paper, we ignore 
these correction terms. 

If we assume degenerate single-particle energies e a = 
0.0, the above Hamiltonian has SO (5) symmetry 01 an d 
the eigenenergy is assigned by the valence nucleon num- 
ber n and isospin T 2j , 



Hlv )sO(5) - 



—Gn 
2 



n - 



n — 6 



(5) 



where ft = J2 a 1S ^ e degeneracy of shell orbits. Note 
that the above equation includes the so-called symmetry 
energy term with coefficient a(A)/A = (G + k°)/2. The 
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FIG. 3: The energy difference between the T — and T = 1 
states in odd-odd N — Z nuclei. The experimental values 
of the differences are denoted by solid diamonds. The open 
squares present the values estimated from the experimental 
odd-even mass differences in Fig. 1 and the T — pairing 
force strength k° . The dashed line is fc° - 10.4yU 1/3 . 



parameters G and k° used above give just the empirical 
symmetry energy formula a(A) = 134.4(1 — 1.52A -1 / 3 ) 
determined by Duflo and Zuker |T^ . 

We next consider energy difference between the lowest 
T = and T = 1 states in odd-odd N = Z nuclei. Odd- 
odd N = Z nuclei with A < 40 have the ground states 
with T = 0, J > except for 34 C1, while the ground states 
of odd-odd N = Z nuclei with 40 < A < 74 are T = 1 
and J = except for 58 Cu. Several authors discussed 
that this degeneracy is attributed to the delicate balance 
between the symmetry energy and pairing correlations, 
and that the energy difference between T = 1 and T = 
states is well reproduced by £r=i — E T=0 = 2a(A)/A — 
2At=i using the value ~ 75 for a(A) and the pairing gap 
A T= i = 12A- 1 / 2 . However, if we substitute the odd- 
even mass difference A„ (Z = even, Z + 1) for At=i, 
the energy difference Et=i — Et=o becomes larger than 
the experimental value. The energy difference can be 
regarded as a measure of competition between the T = 
and T = 1 pairing correlations as seen from the following 
identity, 



T=l 



E- 



T=0 



2(A 



T=0 
pn 



A 



T=l\ 



(6) 



The relationships Aj~° 



k°/2 and A^ =1 



of- 



,( 3 ) 



pn 

fer an alternative relation Et=i ~ Et=o ~ k° — 2AJ, 
for the energy difference except for 30 P and 34 C1. If 
we adopt the parameter k° = 244.5(1 - 1.67 A^ 1 / 3 )/ A 

and the average value of pairing gap 5. 18 A -1 / 3 for A„ , 
we get the dashed curve in Fig. 3, which displays well 
the trend of the experimental values of energy difference 
Et=i —Et=o- Adopting the experimental odd-even mass 
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FIG. 4: The calculated energy differences between the lowest 
T — and T = 1 states in even-even (upper plots) and odd- 
odd (lower plots) N — Z nuclei. The solid diamonds are the 
same as Fig. 3. The open circles denote the energy differ- 
ences obtained by the BCS calculations. The dashed curve is 
2a(A)/A. 

differences for Al 3) and k° = 244.5(1 - 1.67 'A' 1 / 3 )/ A, wc 
obtain the energy difference Et=i — Et=o denoted by the 
open squares. These values nicely reproduce the exper- 
imental values except for 30 P and 34 C1 as shown in Fig. 
3. The disagreements in 30 P and 34 CI are attributed to 
the large deviations of T = pairing gap from the curve 
k°/2 due to the neglect of the shell effects in Fig. 2. 

Moreover, we calculated the T = and T = 1 energy 
differences for odd-odd N = Z nuclei with A > 78, al- 
though there are no experimental data of the energy dif- 
ference. The calculation predicts that 82 Nb has possibly 
the ground state with T = 0, while the other odd-odd 
N = Z nuclei have the T = 1 ground state. We call 
this isospin inversion hereafter. It is well known that 
a similar isospin inversion occurs at 58 Cu. The isospin 



inversion is due to characteristic situation, where the 
Fermi energy lies between large spin and small spin or- 
bits with large energy gap i.e., I/7/2 and 2p 3 / 2 for 58 Cu, 
and lffg/2 and 2p 1 / 2 for 82 Nb. In these cases, the T = 1 
pairing gap is quite small as seen in Fig. 1, and en- 
ergy difference becomes large from the simple relation 
E T=1 - E T=0 w k° - 2Al 3) (Z = even, Z + l). 

Figure 4 shows the calculated energy differences 
Et=i — Et=o in odd-odd and even-even N = Z nuclei. 
The energy differences in the BCS approximations are 
calculated by 2a(A)/A + Abcs f° r even-even N = Z 
nuclei and by k° — 2/S.bcs f° r odd-odd N — Z nuclei 
where a (A) is the empirical symmetry energy coefficient 
and Abcs is the BCS pairing gap. The BCS calculations 
well reproduce the experimental values of energy differ- 
ences, except for odd-odd N = Z nuclei with A < 40. 
The BCS calculations show that the T = and T = 1 
states in 82 Nb are almost degenerate, while the ground 
states of adjacent odd-odd N = Z nuclei have isospin 
T= 1. 

In conclusion, we investigated the T = and T = 1 
pairing correlations in N = Z nuclei. The T = 1 pair- 
ing correlations in N = Z nuclei are extracted from the 
odd-even mass differences of the neighboring odd-mass 
nuclei, which can be fitted by the curve 5.18A~ 1//3 . The 
pn pairing correlations are estimated from the double dif- 
ference of binding energies. The T — 1 pn pairing gap 
is the same as the nn pairing gap. The indicator Aj=° 
presents the magnitude of T = pn pairing correlations. 
The energy differences between the T — and T = 1 
states are well described by the T = 1 and T — pair- 
ing model. In odd-odd N = Z nuclei, the T — 1 pairing 
correlations compete with the T — pairing correlations, 
and the degeneracy of the T = and T = 1 states oc- 
curs. The empirical values and BCS results reproduced 
the energy difference. In particular, our results predict 
that odd-odd N = Z nucleus 82 Nb has the T = ground 
state or the T = and T = 1 states are almost degener- 
ate. The odd-even mass differences for even-even N = Z 
nuclei are extremely larger than those of the neighboring 
even-even N ^ Z nuclei. It would be affected by strong 
pn correlations. Further studies in this direction are in 
progress. 



[1] W. Satula and R. Wyss, Phys. Rev. Lett. 86, 4488(2001); 

87, 052504(2001). 
[2] K. Kaneko and M. Hasegawa, Phys. Rev. C 60, 

024301(1999); Prog. Theor. Phys. 106, 1179(2001). 
[3] W. Satula, J. Dobaczewski, and W. Nazarewicz, Phys. 

Rev. Lett. 81, 3599(1998). 
[4] J. Dobaczewski, P. Magierski, W. Nazarewicz, W. Satula, 

and Z. Szymanski, Phys. Rev. C 63, 024308(2001). 
[5] T. Duguet, P. Bonche, P.-H. Heenen, and J. Meyer, Phys. 

Rev. C 65, 014311(2001). 
[6] J. Janecke, Nucl. Phys. 73, 97(1965). 
[7] N. Zeldes and S. Liran, Phys. Lett. 62B, 12(1976). 



[8] A. O. Macchiavelli, P. Fallon, R. M. Clark, M. Cromaz, 

M. A. Deleplanque, R. M. Diamond, G. J. Lane, I. Y. 

Lee, F. S. Stephens, C. E. Svensson, K. Vetter, and D. 

Ward, Phys. Rev. C 61, 041303(R)(2000); Phys. Lett. 

480B, 1(2000). 
[9] J. Janecke, T. W. O'Donnell, and V. I. Goldanskii, Phys. 

Rev. C 66, 024327(2002). 
[10] S. Frauendorf and J. Sheikh, Nucl. Phys. A645, 

509(1999). 

[11] P. Vogel, Nucl. Phys. A662, 148(2000). 
[12] E. Baldini-Neto, C. L. Lima, and P. Van Isacker, Phys. 
Rev. C 65, 064303(2002). 



■5 



[13] J. Duflo and A. P. Zuker, Phys. Rev. C 52, 23(R)(1995). 
[14] S. Hilaire, J.-F. Berger, M. Girod, W. Satula, and P. 

Schuck, Phys. Lett. 531B, 61(2002). 
[15] J. Janecke and H. Brehrens, Phys. Rev. C 9, 1276(1974). 
[16] A. S. Jensen, P. G. Hansen, and B. Jonson, Nucl. Phys. 

A431, 393(1984). 



[17] M. Hasegawa and K. Kaneko, Phys. Rev. C 61, 
037306(2000); M. Hasegawa, K. Kaneko, and S. Tazaki, 
Nucl. Phys. A688, 765(2000). 

[18] K. T. Hecht, Phys. Rev. 139B, 794(1965); Nucl. Phys. 
A102, 11(1967). 



